Spin current and rectification in one-dimensional electronic systems 
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Spin and charge currents can be generated by an ac voltage through a one-channel quantum wire 
with strong electron interactions in a static uniform magnetic field. In a certain range of low voltages, 
the spin current can grow as a negative power of the voltage bias as the voltage decreases. The spin 
current expressed in units of ft/2 per second can become much larger than the charge current in 
units of the electron charge per second. The system requires neither spin-polarized particle injection 
nor time-dependent magnetic fields. 
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The pioneering paper by Christen and Buttiker [lj 
has stimulated much interest to rectification in quantum 
wires and other mesoscopic systems. Most attention was 
focused on the simplest case of Fermi liquids 0, S| ■ Re- 
cently this research was extended to strongly interacting 
systems where Luttinger liquids are formed 0, M, @, 0] • 
One of the topics of current interest is the rectification ef- 
fect in mesoscopic conductors in a magnetic field 0, S 0| • 
In the presence of a magnetic field, both spin and charge 
currents can be generated. However, only charge currents 
have been studied so far. In this paper we consider the 
generation of a dc spin current by an ac voltage bias in 
a one-channel quantum wire. 

In recent years many approaches to the generation of 
spin currents in quantum wires were put forward. Typ- 
ically both, spin and charge currents are generated and 
the spin current expressed in units of h/2 per second is 
smaller than the electric current in units of e per sec- 
ond (e is the electron charge). Such a situation naturally 
emerges in partially polarized systems since each electron 
carries the charge e and its spin projection on the z-axis is 
±H/2. A proposal how to obtain a spin current exceeding 
the charge current in a quantum wire was published by 
Sharma and Chamon 0, Q who considered a Luttinger 
liquid in the presence of a time-dependent magnetic field 
in a region of the size of an electron wavelength. We show 
that the generation of a dc spin current exceeding the 
charge current is also possible without time-dependent 
magnetic fields on the nanoscalc. The spin current can be 
generated in a spatially asymmetric system in the pres- 
ence of an ac bias. Interestingly, in a certain interval 
of low voltages the dc spin current grows as a negative 
power of the ac voltage when the voltage decreases. 

The rectifying quantum wire is sketched in Fig. [T] It 
consists of a one-dimensional conductor with a scatterer 
in the center of the system at x = 0. The scatterer creates 
an asymmetric potential U{x) ^ U{—x). The size of the 
scatterer au ~ is of the order of the electron wave- 

length. The wire is placed in a uniform magnetic field H. 
The field defines the S z direction of the electron spins. 



At its two ends, the wire is connected to nonmagnetic 
electrodes, labeled by i = 1,2. The left electrode, i = 1, 
is controlled by an ac voltage source, while the right elec- 
trode, i = 2, is kept on ground. The magnetic field H 
breaks the symmetry between the two orientations of the 
electron spin. In a uniform wire this would not result 
in a net spin current since the conductances of the spin- 
up and -down channels would be the same, e 2 /h [llj . 
and the spin currents of the spin-up and -down electrons 
would be opposite. In the presence of a potential barrier 
such cancellation does not occur [12]. In a system with 
strong electron interaction, the spatial asymmetry of the 
wire leads to an asymmetric I —V curve, I(V) ^ J(— V) 
[4|. Thus, an ac voltage bias generates spin and charge 
dc currents, JJ and I r c . In this paper we focus on the 
low-frequency ac bias. We define the rectification cur- 
rent as the dc response to a low-frequency square voltage 
wave of amplitude V: F S {V) = [I S {V) + I s (-V)}/2 and 
I r c {V) = [I C (V) + I c (-V)]/2. The above dc currents ex- 
press via the currents of spin-up and -down electrons: 
I r c = Zf + J[, F s = (V2e)[/f - II]. The spin current 
exceeds the charge current if the signs of IT and 2T are 
opposite. 

The calculation of the rectification currents reduces to 
the calculation of the contributions even in the voltage 
V to the dc I — V curves I S (V) and I C (V). We assume 
that the Coulomb interaction between distant charges is 
screened by the gates. This will allow us to use the stan- 
dard Tomonaga-Luttinger model with short range inter- 
actions [l]| . Electric fields of external charges are also as- 
sumed to be screened. Thus, the applied voltage reveals 
itself only as the difference of the electrochemical poten- 
tials Ei and E% of the particles injected from the left and 
right reservoirs. We assume that one lead is connected to 
the ground so that its electrochemical potential E2 = Ep 
is fixed. The electrochemical potential of the second lead 
Ei = Ep + eV is controlled by the voltage source. Since 
the Tomonaga-Luttinger model captures only low-energy 
physics, we assume that eV -C Ep, where Ep is of the 
order of the bandwidth. 
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Rectification occurs due to backscattering off the 
asymmetric potential U{x). We will assume that the 
asymmetric potential is weak, U{x) -C Ep, This will 
enable us to use perturbation theory. As shown in Ref. 

the leading contribution to the rectification effect 
emerges in the third order in U (x). This fact allows for a 
simple explanation of our main prediction that the spin 
rectification current can significantly exceed the charge 
current. To illustrate the principle, let us consider the 
following toy problem: There is no uniform magnetic 
field H and no asymmetric potential U{x). Instead, 
both right<-deft and spin-up<->spin-down symmetries are 
broken by a weak coordinate-dependent magnetic field 
B z (x) B z (—x), which is localized in a small region of 
size ~ 1/kp (we do not include the components B XiV in 
the toy model). Let us also assume that the spin- up and 
-down electrons do not interact with the electrons of the 
opposite spin. Then the system can be described as the 
combination of two spin-polarized one-channel wires with 
opposite spin-dependent potentials ±/j,B z (x), where /z is 
the electron magnetic moment. According to Ref. [4{ 
an ac bias generates a rectification current in each of 
those two systems and the currents are proportional to 
the cubes of the potentials (±{iB z ) 3 . Thus, II = —IT 



nian would further be diagonalized by the combinations 
4> c .s oc (f>i ±</>j., and similarly for H c , s , expressing the spin 
and charge separation. This is here no longer the case. 
If we focus on the cj) fields only (as II will not appear 
in the operators describing backscattering off U(x)), the 
fields diagonalizing the Hamiltonian, 4> c ,s, have a more 
complicated linear relation to <£t,j., which we can write 
as 



a] 



P] 
-0] 



(2) 



and which corresponds to the matrix A of [12j | . The nor- 
malization has been chosen such that the propagator of 
the <j> fields with respect to the Hamiltonian (fTJ) evaluates 
to (<f)c,s{ti)(t>c,s(t2)) = — 21n(z(ti— t 2 ) + S), where 5 > is 
an infinitesimal quantity. For non-interacting electrons 
without a magnetic field, g c = g s = 1/2. g c < 1/2 
(> 1/2) for repulsive (attractive) interactions. In the 
absence of the magnetic field, terms in (fTJ) in the form 
of exp(±2iy / 570 s ) may become relevant and open a spin 
gap for g s < 1/2. In our model they can be neglected 
since they are suppressed by the rapidly oscillating fac- 
tors exp(±2i[fci?| — kpi\x). It is convenient to model the 



Hence, no net charge current V = IT + IT is generated in leads as the re S ions near the rigMand left ends of the 



the leading order. At the same time, there is a nonzero 
spin current in the third order in B z . 

Let us now calculate the currents in the presence of the 
asymmetric potential U(x) and the field H. We assume 
that the magnetic field couples only to the electron spin 
and neglect the correction — eA/c to the momentum in 
the electron kinetic energy. Indeed, for a uniform field 
one can choose A ~ y, where the y-axis is orthogonal to 
the w ire, and y is small inside a narrow wire. As shown in 
Ref. [l2j, such a system allows a formulation within the 
bosonization language and, in the absence of the asym- 
metric potential, can be described by a quadratic bosonic 
Hamiltonian 



wire without electron interaction [11 1 

Backscattering off the impurity potential U(x) is de- 
scribed by the following contribution to the Hamiltonian 
H = H Q + H> fli: 



(3) 



H 



E 



E 

vy=L,Ra,a' = ],l 



dx{d x (j) va )'H va y a ' {d x 4> V 'a>), 



(1) 

where a is the spin projection and v — R, L labels the 
left and right moving electrons, which are related to the 



boson fields 



as 



± for v — R, L. The operators are the Klein factors 
adding a particle of type (v, a) to the system, and kp va j-K 
is the density of (v, a) particles in the system. The densi- 
ties of the spin-up and -down electrons are different since 
the system is polarized by the external magnetic field. 
The 4x4 matrix H describes the electron-electron inter- 
actions. In the absence of spin-orbit interactions, L «-> R 
parity is conserved and we can introduce the quantities 
4><r = 'pLa+'t'Ro- and = 4>l<j—4>r.<t such that the Hamil- 
tonian decouples into two terms depending on <f> a and n CT 
only. In the absence of the external field, this Hamilto- 



where the fields are evaluated at position x — and 
U (nj, ?i|) = U*{— n-f, — njj since the Hamiltonian is Her- 
mitian. The fields n do not enter the above equation 
due to the conservation of the electric charge and the 
z-projection of the spin. The Klein factors are not writ- 
ten because they drop out in the perturbative expansion. 
U(n^ , n± ) are the amplitudes of backscattering of n-\ spin- 
up and rt| spin-down particles with n a > for L — > R 
and Ti a <C for R — ► L scattering. U{n^ni) can be 
estimated as £/(n T ,raj,) - / dxU(x)e m ^ k n x+m t 2k Fi x ~ 
U/kp, where U is the maximum of U(x) (cf. (HI). In 
the case of a symmetric potential U(x) = U(—x) the 
coefficients U(n^,n^) are real. 

The spin and charge current can be expressed as 



I s ,c — L sc 



nl _ r2 



(4) 



where L l s c and R\ c denote the current of the left- and 
right-movers near electrode i, respectively. For a clean 
system (U(x) = 0), the currents Rl c — R^ c , L\ c = 
L 2 S C and I c = 2e 2 V/h, I s = With backscattering 
off U(x), particles are transferred between L and R in 
the wire, and hence R 2 = + dSp/dt, R 2 C = R\ + 
dQji/dt, where Qr and Sp denote the total charge and 
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the z-projection of the spin of the right-moving electrons 
14J. The currents L 2 C and R\ c are determined by the 
leads (i.e. the regions without electron interaction in our 
model [llj]) and remain the same as in the absence of the 
asymmetric potential. Thus, the spin and charge current 
can be represented as I c = 2e 2 V/h + I bs and L = I bs , 
where the backscattering current operators are [ISO 



4 =dQ R /dt = ie[H,Q R ]/h 



= ~jf E ( n T + n l) U ( n h n l) e 



in T T (O)+m i i (O) 



(5) 



I b s s = dS R /dt = - 5 E (™T - n i )U(n h n l )e m ^ +m ^K 

(6) 



The calculation of the rectification currents reduces to 
the calculation of the currents ([5]), (O at two opposite 
values of the dc voltage. 

To find the backscattered current we use the Keldysh 
technique (la ]. We assume that at t = — oo there is 
no backscattering in the Hamiltonian (U(x) — 0), and 
then the backscattering is gradually turned on. Thus, 
at t = — oo, the numbers Np and N R of the left- and 
right-moving electrons conserve separately: The system 
can be described by a partition function with two chem- 
ical potentials E\ — Ep + eV and E2 = Ep conjugated 
with the particle numbers N R and Np. This initial state 
determines the bare Keldysh Green functions. 

We will consider only the zero temperature limit. It 
is convenient to switch (l3 | to the interaction represen- 
tation Ho — ► Ho — EiN R — E2N1,. This transformation 
induces a time dependence in the electron creation and 
annihilation operators. As the result each exponent in 
Eq. ((3]) is multiplied by exp(ieVt[n^ + n\\/K). 

In the Keldysh formulation [l5| the backscattering cur- 
rents ([5]), ([6]) are evaluated as 



C = (0|5(-oo,0)/^5(0,-oo)|0), 



(7) 



where |0) is the ground state for the Hamiltonian Hq, Eq. 
(TTJ), and S(t,t') the evolution operator for H' from t' to 
t in the interaction representation with respect to Hq. 
The result of this calculation depends on the elements of 
the matrix ([2|) , which describe the low-energy degrees of 
freedom and depend on the microscopic details. Several 
regimes are possible id] at different values of the param- 
eters <7 S > 0, g c > 0, a and f3. In this paper we focus on 
one particular regime: 



1 > g c > g s > g c - g s > a, (3. 



(8) 



We will see that in such a regime the spin current can be 
much greater than the charge current. We will also as- 
sume that [g c ot + g s (3] > but our results do not depend 
significantly on this assumption. The behavior which we 



find in the regime (8) persists for g c ,g s ,a and j3 smaller 
than some constants of the order of unity. In this pa- 
per we focus on the limit in which the current can be 
calculated analytically. A numerical investigation of the 
boundary of the region where the spin current exceeds 
the charge current is beyond the scope of the present 
paper. 

The currents ([7|) can be estimated using a renormal- 
ization group procedure [l3j |. As we change the energy 
scale E, the backscattering amplitudes U(n^, n^) scale as 
U(n : m; E) ~ E zl ^ n,m \ where the scaling dimensions are 

z{n,m) = n 2 [g c (l + a) 2 + g s (l + (3) 2 ] + m 2 [g c {\ - af 
+g s (l - I3f] + 2nm[g c (l - a 2 ) - g s (l - p 2 )} - 1.(9) 



The renormalization group stops on the scale of 
the order E ~ eV. At this scale the backscatter- 
ing current can be represented as I bs s = Vr c 
where the effective reflection coefficient r C)S (V) 
is given by the sum of contributions of the form 
(const) [/(m, mi)U(n 2 , m 2 ) • . . U(n p , m p )V^=i z (™*> m *). 
Such a perturbative expansion can be used as long 
as U is small. In the regime (JSJ) it is sufficient to 
require that U(n,m) < (eV) 1+s /E s F , where S > g c . 
The contributions to the backscattering current of the 
nth order in U scale as I^) ~ jj n y 1 ~ n +o(gc) _ Hence, 
they exceed the contributions of the (n + l)th order 
for U(n,m) < (eV) 1+s /E F . The leading non-zero 
contribution to the backscattered current, therefore, 
emerges in the second order. The rectification current, 
however, is dominated by a third order contribution for 
V 1+s > U > V 2 /E F . 

Indeed, the second order contributions to the charge 
current were computed in Ref. [131 ]. The spin current 
can be found in exactly the same way. The result is 

lP s bs (V) - ^(const)|J7(n,m)| 2 |F| 22 ("' m ) +1 signK 

(10) 

If the U(n, m) were independent of the voltage, the above 
current would be an odd function of the bias and hence 
would not contribute to the rectification current. The 
backscattering amplitudes depend [l3| on the charge den- 
sities kp va though, which in turn depend on the volt- 
age in our model [4[. The voltage-dependent corrections 
to the amplitudes are linear in the voltage at low bias. 
Hence, the second order contributions to the rectification 
currents scale as u 2 \V\ 2z ^ n,m " l+2 . The additional factor of 
V makes the second order contribution smaller than the 
leading third order contribution (|13jl at U ^> V 2 /Ep. 
Note that the second order contribution to the rectifi- 
cation current is nonzero even for a symmetric poten- 
tial U(x). The leading third order contribution emerges 
solely due to the asymmetry of the scatterer. 

It is easy to find the most relevant backscattering op- 
erators in the renormalization group sense using Eqs. ([8|) 
and ©• The most relevant operator is U(1,0), the sec- 
ond most relevant {7(0, —1), and the third most relevant 
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U(— 1, 1). The leading non-zero third order contributions 
to the spin and charge currents come from the product 
of the above three operators in the Keldysh perturbation 
theory. This leads to 



r. 



jy3y 9c (2+6a 2 )+ Ss (6+2 / 3 2 )-2 



(11) 



This contribution dominates at (eV) 2 /Ep <C U <C (eV). 
Interestingly, the current (jTTJ) grows as the voltage de- 
creases. 

Does the current (jlip contribute to the rectification 
current? In general, (jlip is the sum of odd and even func- 
tions of the voltage and only the even part is important 
for us. One might naively expect that such a contribution 
has the same order of magnitude for the spin and charge 
currents. A direct calculation shows, however, that this 
is not the case and the spin rectification current is much 
greater than the charge rectification current. 

In order to calculate the prefactors in the 
right hand side of Eq. (fTTj) one has to em- 
ploy the Keldysh formalism. The third order 
Keldysh contribution reduces to the integral of 
P(ti,t 2 ,t 3 ) = (Tcexp^ftiJ + ieVti/fiJeBtpH^^)- 
ieVt 2 /h)exp(i{-<f> 1 (t 3 ) + </>i(f 3 )D) over (t x - i 3 ) and 
(*2 — is), where T c denotes time ordering along the 
Keldysh contour — oo — ► — ► — oo and the angular 
brackets denote the average with respect to the ground 
state of the non- interacting Hamiltonian fTj). For the 
purpose of comparing the prefactors in Eq. (jTTJ) for 
and Ig S in the region of parameters (j8j), it is sufficient 
to perform the calculation in the limit 1 » j c = j s . 
a = = 0, which corresponds to the boundary of the 
region ((SJ. A small change of the parameters g c ,g s ,a 
and f3 will not significantly affect the prefactors (as well 
as the exponent in Eq. (fTTj) ). In the above limit the 
calculation considerably simplifies and can be performed 
analytically. The correlation function P(ti,t2,ts) fac- 
torizes into the product of a function of (t\ — £3) and a 
function of [to, —t$). One finds 



h 3 



sign(eF) 



eVr r 



r 2 (l-4.g)sin(47r 5 ) 



x [1 - cos(47rff)]Re[Z7(l, 0)Z7(— 1, 1)U(0, -1)], 



(12) 



jbs 



At 2 

_ ' C 

h 2 



eVr c 



8g~2 



r 2 (i-4. 9 ) 



x sin 2 (47r ff )Im[;7(l, 0)[/(-l, l)C/(0, -1)], (13) 

where g = g c = g s and r c ~ h/ Ep is the ultraviolet cutoff 
time. The above results apply at g < 1/4, i.e., when the 
Keldysh integral converges. The charge current (jT2|) is an 
odd function of the voltage and hence does not contribute 
to the rectification effect. The spin current (fT3|) is an even 
function and hence determines the spin rectification cur- 
rent. It is non-zero if Im[U(l, 0)f7(-l, 1)17(0, -1)] ^ 0, 




quantum wire 



U(x) 



au 




FIG. 1: Sketch of the one-dimensional conductor connected 
to two electrodes on both ends. Currents are driven through 
a voltage bias V that is applied on the left electrode while the 
right electrode is kept on ground. The system is magnetized 
by the field H. Electrons are backscattered off the asymmetric 
potential U(x). U(x) 7^ in the region of size au ~ I/Uf- 



which is satisfied for asymmetric potentials. Thus, we ex- 
pect that in the region ([5]), the spin rectification current 
exceeds the charge rectification current in an appropriate 
interval of low voltages \fU Ep 3> eV 3> U. 

The voltage dependence of the spin rectification cur- 
rent is illustrated in Fig. [2j The expression (fTTj) describes 
the current in the voltage interval V** > V > V* . In this 
interval the current increases as the voltage decreases. 
At lower voltages the perturbation theory breaks down. 
The current must decrease as the voltage decreases be- 
low V* and eventually reach at V = 0. At higher volt- 
ages, E F 3> eV ^> eV** , the second order rectification 
current (jlOj) dominates. The leading second order con- 
tribution ij ~ \U(l,0)\ 2 V 2z ^ +2 grows as the voltage 
increases. The charge rectification current has the same 
order of magnitude as the spin current. The Tomonaga- 
Luttinger model cannot be used for the highest voltage 
region eV ~ Ep. 

In conclusion, we have shown that rectification in mag- 
netized quantum wires can lead to a spin current that 
largely exceeds the charge current. The effect is solely 
due to the properties of the wire and does not require spin 
polarized injection as from magnetic electrodes. The cur- 
rents are driven by the voltage source only. In an interval 
of low voltages the spin current grows as the voltage de- 
creases. 
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